Self Consistent Simulation of Quantum Transport and Magnetization Dynamics in 

Spin- Torque Based Devices 
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We present a self consistent solution of quantum transport, using the Non Equilibrium Green's 
Function method, and magnetization dynamics, using the Landau-Lifshitz-Gilbert formulation. We 
have applied this model to study 'spin torque' induced magnetic switching in a device where the 
transport is ballistic and the free magnetic layer is sandwiched between two anti-parallel (AP) 
ferromagnetic contacts. We predict hysteretic current-voltage characteristics, at room temperature, 
with a sharp transition between the bistable states that can be used as a non-volatile memory. 
We further show that this AP penta layer device may allow significant reduction in the switching 
current, thus facilitating integration of nanomagnets with electronic devices. 

PACS numbers: 72.25.Dc 



Successful integration of nanomagnets with electronic 
devices may enable the first generation of practical spin- 
tronic devices, which have so far been elusive due to strin- 
gent requirements such as low temperature and high mag- 
netic field. It was predicted by Slonczewskii and Berges^ 
that magnetization of a nano magnet may be flipped by a 
spin polarized current through the so-called 'spin torque' 
effect and this was later demonstrated experimentallji^. 
However, the early spin-torque systems were metal based 
that allowed only a small change in the magnetoresis- 
tance. In addition, metallic channels are difficult to in- 
tegrate with CMOS technology. Recently a number of 
experiments have demonstrated current-induced magne- 
tization switching in MgO based Tunneling Magneto Re- 
sistance (TMR) devices at (i) room temperature (ii) with 
a TMR ratio of more than 100% and (iii) without any 
external magnetic field^. Encouraged by these experi- 
mental results, here we explore theoretically a memory 
device based on current induced magnetization switching 
in the quantum transport regime. 

The device under consideration is shown in Fig. ^ 
It consists of five layers. The two outer layers are 'hard 
magnets' which act as spin polarized contacts. There is a 
soft magnetic layer inside the channel whose magnetiza- 
tion is affected by the current flow through the so-called 
'spin torque' effect. The channel can be a semiconductor— 
or a tunneling oxida&. Note that the contacts are ar- 
ranged in an anti-parallel (AP) conflguration. We have 
recently showed that in this configuration, the torque ex- 
erted by the injected electrons on a the nearby spin array 
(in this case the soft magnet) is maximum^. A similar 
prediction was also made by Berger— based on expan- 
sion/contraction of the Fermi surface. The possibility of 
an enhanced torque and therefore a lower switching cur- 
rent is our motivation for the penta-layer configuration 
instead of the conventional tri-layer geometry. 

In Fig. 1, the soft-magnet changes the transport 
through its interaction with the channel electrons, which 
in turn exert a torque on the magnet and try to rotate 
it from its equilibrium state. In this paper, we present 



a self-consistent solution of both these processes: the 
transport of channel electrons (through NEGF) and the 
magnetization dynamics of the free layer (through LLG 
equations) (see Fig. H^b)). Our calculations show clear 
hysteretic I-V suggesting possible use as a memory. Fur- 
thermore, we show that a penta-layer device with AP 
contact as shown in Fig. 1(a) should exhibit a significant 
reduction in the switching current. 

Transport.- Unlike the conventional metallic spin- 
torque systems, where the transport is predominantly 
diffusive, the transport in semiconductors or tunneling 
oxides is ballistic or quasi-ballistic. This necessitates a 
quantum description of the transport. We use the Non 
Equilibrium Green's Function (NEGF) method to treat 
the transport rigorously. The interaction between chan- 
nel electrons and the ferromagnet is mediated through ex- 
change and it is described by Hj{r) = J{r—Rj)(j.Sj, 
where, r and Rj are the spatial coordinates and a and 
Sj arc the spin operators for the channel electron and 
j-th spin in the soft-magnet. J(f — Rj) is the interaction 
constant between the channel electron and the j-th spin 
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FIG. 1: (a) A schematic showing the penta layer device. The 
free ferromagnetic layer is embedded inside the channel which 
is sandwiched between two 'hard' ferromagnetic contacts, (b) 
A schematic showing the self-consistent nature of the trans- 
port problem. The magnetization dynamics and transport are 
mutually dependent on one another. 
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FIG. 2: Non self consistent (with magnetization dynamics) pj^ 3. ^j^^ hysteretic I-V originating from a self-consistent 

I-V characteristics of the the proposed device (a) With the solution of transport and LLC. At a certain bias, the current 

soft magnet initially at ^ = ^ position. THe current is larger ^^^^^^ produced by the conduction electrons is strong enough 

for positive bias, (b) With the soft magnet initially at = ^g^^riet. These transition points are indicated in 

position. The current is larger for negative bias. . v fieure 



in the magnet. This interaction is taken into account 
through self energy (Es), which is a function of the mag- 
netization (m), using the so-called self-consistent Born 
approximationifl. In this method, the spin current flow- 
ing into the soft magnet is given by 

[Lp^n] = J dE^i [Tr {GSr - J^TG^ - S,G" + G"I]t}] 

(1) 

where the trace is taken only over the space coordi- 
nates. Then [/spin] is a 2 x 2 matrix in the spin space. 
Here, G denotes the Green's function. The torque ex- 
erted on the magnet is calculated from [/spin] by writ- 
ing Ti = Trace {Si [/sp,„]}, where i = {x,y,z}. The to- 
tal current, which is found from a similar expression as 
Eq. n with the self energy now replaced by the total 
self-energy Sii, is shown in Fig|2Ia),(b) for two differ- 
ent configuration of the magnetization. The nonlinear 
nature of the I-V can be intuitively understood by rec- 
ognizing that the exchange interaction is minimum if the 
injected electrons and the magnetization has the same 
spin orientation^. 

For the calculations, the Hamiltonian was written in 
the effective-mass approximation where the hopping pa- 
rameter t = h^/{2m*a?), TO* = O.TTOjfi^ denoting the 
effective mass and a being the lattice spacing. The in- 
teraction constant J is assumed to be 0.01 eV— . A sam- 
ple set of parameters is = 2eV ; barrier-height=1.2 
e\f—] barrier- width=l nm and exchange splitting=1.2 
eV—. However, the barrier height, width and exchange 
splitting were artificially varied to get desired injection 
efficiency and TMR value. We also modified the perfect- 
contact self energy to read T, = —t exp{ika){t ^ t and 
k is the momentum) to simulate the reflective nature of 
the contact (for details seai^). For plots 2-4, an injection 
efficiency of 70% was assumed. 

Magnetization Dynamics. -The magnetization dynam- 
ics is simulated using the LLG equation: 
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FIG. 4: Response to a Read-Write-Read pulse. The write 
pulse switches the magnet from 9 — n to 9 — 0. The corre- 
sponding change in the current can be clearly seen during the 
write pulse. 



Here, m is the magnetization of the soft magnet, 7 = 
17.6 MHz/Oe is the gyro magnetic ratio and a is the 
Gilbert damping parameter. The H^jf = H^xt + 
{2Ku2/Ms)mzZ — {2Kup/Ms)mxX, where Hext is the 
externally applied magnetic field, Ms is the saturation 
magnetization and Ku2 and Kup are the uniaxial and 
in-plane anisotropy constants respectively. The conven- 
tional LLG equation has to be solved with the current 
torque (Ti) that works as an additional source term. 

Self-Consistency. -Fig. [21shows the situations when the 
transport and magnetization dynamics is independent of 
each other. This will change when Eqs. n]and|2]are solved 
self-consistently. If we start from 9 = it position, I-V 
curve follows the trend shown in FiglJl^a). However, once 
the torque exceeds the critical field (discussed later), the 
magnet switches abruptly. As a result I-V characteristics 
now follows that shown in FigI2Ib). This results in the 
hysteretic I-V shown in Fig. O 

Fig^ shows current flow in the device in response to 
Read-Write-Read pulse sequence. Here, we have used 
read pulse of 0.5 volt and write pulse of -|-1 volt. The soft 
magnet is initially in the 9 = tt position. The Write pulse 
switches it to 6* = 0. Note the change in the current level 
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in response to the Read pulse before and after applying 
the Write pulse. 

Discussion. -A question may be raised regarding the 
asymmetric I-V of Fig. |3 which is not expected if one 
thinks about the device in Fig. Q^a) as a series com- 
bination of two devices, one anti-parallel (AP) and one 
parallel (P). The device, however, is different from a mere 
series combination since the contact in the middle works 
as a mixing element for up and down spin electrons. The 
difference will be clear if one assumes 100% injection effi- 
ciency. No current is expected to flow through the series 
combination of an AP and a P device. However, in our 
device, a current can still flow because the contact in the 
middle mixes the up and down spin channels. This 'ex- 
tra' current originating from 'channel-mixing' gives the 
observed asymmetry in Fig. |21 

Since electronic time constants are typically in the sub 
picosecond regime which is much faster than the magneti- 
zation dynamics (typically of the order of nano seconds), 
we have assumed that, for electronic transport, the mag- 
netization dynamics is a quasi-static processi^. 

The switching is obtained by the torque compo- 
nent which is transverse to the magnetization of the 
soft-magnet. From Eq. [21 considering average rate 
of change of energy, it can be shown that the mag- 
nitude of the torque required to induce switching is 
{H^^t + Hk + where Hu = 2Ku2/Ms and 

Hp = 2Kup/Ms = AttMs. This then translates into a 
critical spin current magnitude of 

Isp^n = ^a{M,V){H,,t + Hk + 2nA'Q. (3) 

Here, V is the volume of the free magnetic layer. De- 
pending on the magnitudes of a,Mg,H}^ and thickness, d, 
of the magnet, the spin current density to achieve switch- 
ing varies from 10^ — 10^ A/cm^ (e.g for Co, using typ- 
ical values a ~ 0.01,7Jfe ~ 100 Oe and = 1.5 x 10^ 
emu/cm'^i^ and d = 2 nm, the spin current density re- 
quired is roughly 10^ A/cm^). Note that this require- 
ment on spin current is completely determined by the 
magnetic properties of the free layer. The actual current 
density is typically another factor of 10-100 larger due to 
the additional coherent component of the current which 
does not require any spin-flip. Hence an important met- 
ric for critical current requirement is r = I coherent / 1 spin , 
which should be as small as possible. Intuitively, with 
AP contacts, the coherent current Icoherent cc t^2a/?, 
where t is the hopping matrix element, a is the ma- 
jority(minority) density of states for the injecting con- 
tact and (3 is the minority(majority) density of states 
of the drain contact. Similarly the spin- flip current 
Isf DC [q:^(1 — Pa) — P'^Pa] , whcrc Pa is the proba- 
bility of a spin in the free layer to be in state o^. It 
follows that 

r , +2 1 _ p2 

-^coherent i ^ ^ c / 1\ 

= -^^""^ = ^ + (i - Pa)il + P^) ' 

where, Pc = {a — (3)/{a + (3) indicates the degree of 
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FIG. 5: Variation of the ratio of Icoherent/ 1 spin for a 3 layer 
and AP penta layer spin-torque device showing the possible 
reduction of switching current for the AP penta layer device 
compared to the 3-layer device. 



contact polarization. These approximate analytical ex- 
pressions (Eqs. 21) agree quite well with detailed NEGF 
calculations described above. The Icoherent ^tnd I spin can 
be found respectively from the symmetric and asymmet- 
ric portion of the non-linear I-V shown in Fig. 2. Fig. 

5 shows the variation of g = rtn-layer /tAP-penta-layer 

with Pc- The plot shows that g » 1 for reason- 
able values of Pc, indicating a lower switching current 
for the penta-layer device. Recent experiments on AP 
penta-laycr devicesiSiiSiSii have shown similar reduction 
of switching current compared to 3-layer devices. These 
experiments seem to follow the general trends of Fig. 5 
as the reduction factor is seen to increase with increasing 
TMR (see Fig. 4 oiiS). A detailed study of the depen- 
dence of the reduction factor on material parameters is 
beyond the scope of this paper. 

The sharp transition between HIGH and LOW states 
in Fig. 3 arises from the bistable nature of the solutions 
to the LLG equation in the absence of any external field 
perpendicular to the easy axis. The intrinsic speed de- 
pends on cij = where B can be roughly estimated 
as i? '--^ hT / {2^b)- a higher speed will require higher 
current density. 

In conclusion, we have shown a scheme for calculat- 
ing the 'spin-current' and the corresponding torque di- 
rectly from transport parameters within the framework 
of NEGF formalism. A non-linear I-V is predicted for 
AP penta-layer devices. Experimental observation of 
this non-linearity (which can also be detected as steps 
or peaks in respectively the first and second derivative of 
the I-V^) would provide strong confirmation of our ap- 
proach. We have further coupled the transport formal- 
ism with the phenomenological magnetization dynamics 
(LLG equation) . Our self consistent simulation of NEGF- 
LLG equations show clear hysteretic switching behavior, 
which is a direct consequence of the non-linearity de- 
scribed above. Finally, we have shown that the switch- 
ing current for AP penta layer devices can be significantly 
lower than the conventional 3 layer devices. 
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